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I. INTRODUCTION 

Binary black holes are among the most promising 
sources of gravitational radiation for the Laser Inter- 
ferometer Gravitational wave Observatory (LIGO) and 
other gravitational wave interferometers. LIGO has re- 
cently reached its design sensitivity, making the detection 
of astrophysical sources of gravitational radiation a dis- 
tinct possibility. To aid both in the identification of such 
a signal and in its interpretation, theoretical templates of 
potential gravitational wave forms are urgently needed. 

Numerical relativity is the most promising tool for 
modeling the coalescence and merger of binary black 
holes. Typically, numerical relativity calculations com- 
pute a solution to Einstein's equations in two steps. In 
the first step, the constraint equations of general relativ- 
ity are solved to construct initial data, describing a snap- 
shot of the gravitational fields at one instant of time, and 
in the second step these initial data are evolved forward 
in time by solving the evolution equations. 

Different approaches have been used to construct bi- 
nary black hole initial data (see, e.g., the reviews P, I3 
and references therein). Most earlier calculations em- 
ployed the so-called transverse-traceless decomposition 
of the constraint equations (e.g. H, IE Hli see also 
below), while most more recent calculations solve the 
constraint equations in the conformal thin-sandwich for- 
mahsm (e.g. 0, H H [13; see also [H El El H for 
alternative ways of constructing binary black hole initial 
data.) There is general consensus that the latter formal- 
ism is better suited for the construction of quasiequilib- 
rium data (but see (15t | for a very promising alternative 
approach), even though, at least in terms of global quan- 
tities, both formalisms lead to very similar results for 
configurations outside the innermost stable circular orbit 
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(see, e.g., HEl)- 

Dynamical evolutions of binary black holes have long 
suffered from numerical instabilities. The past year, how- 
ever, has seen dramatic progress, and at this point inde- 
pendent codes using very different approaches and tech- 
niques can reliably model the coalescence, merger and 
ring-down of binary black holes (17, 18, ISjJfl,,^]. Some 
of these calculations 0, 0, |^ treat the black hole sin- 
gularity with the so-called puncture method (see also 
|23L |2J|). Given that puncture initial data for binary 
black holes are most easilyconstructed in the transverse- 
traceless formalism (see Q and compare [H H^), this 
development has renewed some interest in transverse- 
traceless initial data. 

In this paper we construct black hole binary initial 
data for arbitrary mass ratios perturbatively. We adopt 
the transverse-traceless decomposition together with the 
puncture approach and treat both the effect of each black 
hole's boost as well the effect of the companion as a per- 
turbation of a spherically symmetric Schwarzschild black 
hole. Both effects are axisymmetric -albeit with differ- 
ent axes of symmetry- allowing us to find simple ana- 
lytic expressions for the gravitational fields, the location 
of the apparent horizon, the irreducible mass, and the 
total energy. To leading order we can construct binaries 
by simply adding the individual corrections, resulting in 
analytic, perturbative black hole initial data. The solu- 
tion is accurate up to order (P/M)^ and (M/s), where P 
is the black hole's momentum, M a measure of the black 
hole's mass, and s the binary separation, and becomes ex- 
act in the limit of infinite separation. We find Newtonian 
expressions for the energy. While this is not surprising, it 
is certainly reassuring. Recovering the Newtonian limit 
is also non-trivial, since it requires taking into account 
the effect of both the boost of the black holes, and the 
distortion due to the companion black hole. 

The numerical construction of initial data for bi- 
nary black holes in quasi-circular orbit requires signif- 
icant computational resources, since it involves step- 
ping through nested iterations of elliptic solves, root- 
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finding and sometimes even minimizations along certain 
sequences (see 0, 0, ^ and Section liii Dl below) . Our 
perturbative framework allows for a completely analyt- 
ical treatment, making the construction of approximate 
binary black hole initial data for arbitrary mass ratios re- 
markably simple. Our results may therefore be of interest 
as initial data for dynamical simulations, especially for 
large binary separations where the errors are relatively 
small. Perhaps more importantly our results provide an- 
alytical insight into the structure of black hole binaries. 
Several of our intermediate results have been derived pre- 
viously (e.g. |23,|2l|2i|33), but to the best of our knowl- 
edge they have never been combined to construct binary 
black hole initial data. 

Our paper is organized as follows. In Section ^ we 
briefly review the initial value problem of general rela- 
tivity and the transverse-traceless decomposition of the 
constraint equations in particular. We present perturba- 
tive solutions of boosted black holes, of black holes with 
companions, and finally of binaries in quasicircular orbit 
in Section ITTTI In Section ITvl we outline how the results 
of the previous Sections can be used to construct ap- 
proximate binary black hole initial data for black holes 
of arbitrary mass ratio in quasicircular orbit. We con- 
clude with a discussion and summary in Section We 
also provide several Appendices that contain all deriva- 
tions and details omitted in the main body of the text. 
Throughout this paper we adopt geometrical units in 
which G = c = 1 . 



II. THE INITIAL VALUE PROBLEM 



Under a 3-1-1 decomposition l31|,|32|, Einstein's equa- 
tions of general relativity split into a set of constraint 
equations -the Hamiltonian constraint and the Momen- 
tum constraint- and a set of evolution equations. Con- 
structing a set of initial data requires specifying a spa- 
tial metric 7ij and an extrinsic curvature Kij on a spa- 
tial hypersurface E that satisfy the constraint equations. 
Such solutions are usually constructed with the confor- 
mal method, whose most general form is given in [33 . l34j . 
Here, we investigate a special case of the general formal- 
ism, which was known earlier, and which is amenable to 
analytical treatment. We introduce a conformal trans- 
formation of the spatial metric 



(1) 



where ip is the conformal factor and -jij the conformally 
related background metric. We also split the extrinsic 
curvature Kij into its trace K and a conformally rescaled 
trace-free part Aij according to 
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(2) 



In terms of these variables the Hamiltonian constraint 
becomes 



8VV -ijR- -iJj^K'^ + tp-'^AijA'^ = 0, (3) 



and the momentum constraint 



(4) 



Here = j^^DiDj is the Laplacian, Di the covariant 
derivative, and R the Ricci scalar associated with the 
background metric 7^ , and we have also assumed vac- 
cuum. 

Both the conformal background metric 7^ and the 
trace of the extrinsic curvature K remain freely speci- 
fiable. We choose conformal flatness 7^ — fij, where fij 
is the flat metric, and maximal slicing K = 0, in which 
case the constraint equations reduce to 
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and 



DjA'^ = 0. 



(5) 



(6) 



Quite remarkably, the momentum constraint ((HJ becomes 
linear and decouples from the Hamiltonian constraint ^ . 
An analytical "Bowen-York" solution to the momentum 
constraint, describing a black hole at coordinate location 
C with linear momentum P is given by [s^ |3^ l37| 



^CP 



2rl 



PkTl 



where re = || 
the center of the black hole and nl 



(7) 

C'll is the coordinate distance from 
{x' - /re is 
the unit vector pointing from that center to coordinate 
location x. Given the linearity of the momentum con- 
straint we can construct solutions for two boosted black 
holes by simple superposition. 



^ — ^CiPi 



^C2P2 



(8) 



Equations {7)) or © are now substituted back into the 
Hamiltonian constraint Eq. JSJ, which is then solved for 
the conformal factor. This step requires boundary con- 
ditions, which must enforce the existence of black holes 
in the constructed initial data. In this paper, we use the 
puncture method (see also [23I to accomplish 
this. In this approach, an appropriate (singular) piece is 
split off the conformal factor analytically, and Eq. (|3J) is 
rewritten as an elliptic equation for the remainder, which 
is continuous and finite throughout. Typically, this latter 
equation is solved numerically. We will instead construct 
an approximate but analytical solution perturbatively, as 
described in the following Section. 
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III. PERTURBATIVE SOLUTIONS 



where 



We construct perturbative binary black hole solutions 
by separately considering the effects of the boost and 
the companion on each black hole, and then combining 
the results. In Section IlII Al we first consider an isolated 
black hole of bare mass Ai with boost P and construct 
a solution that is accurate to second order in ep, where 



P 



(9) 



In Section IIII Bl we then consider a static black hole of 
bare mass A4i with a companion of bare mass A4 2 s^t sl 
coordinate distance s and find a solution that is accurate 
to first order in ec, where 



M 



(10) 



In order to avoid unnecessarily complicated notation we 
do not distinguish between A4i and M2 in ep and e^, 
and instead simply point out that these expressions do 
depend on the mass ratio. 

Finally, in Section flll ("I we combine results to find per- 
turbative binary black hole solutions. For systems in 
equilibrium the Virial theorem implies that we must have 



4, 



so that the orders of expansion for our treatment of the 
boost and the companion are consistent with each other. 



A. A Single Boosted Black Hole 

Consider a static black hole with bare mass A4 at co- 
ordinate location C. The asymptotically flat solution to 
the Hamiltonian constraint © is then the well-known 
expression 
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M 
2rc' 



(12) 



describing a Schwarzschild black hole in isotropic coordi- 
nates. 

To generalize this solution for a boosted black hole 
with momentum P we adopt the puncture method and 
split the conformal factor into two terms. 



1p = h Up. 

a 



Here 



1 

- = 1 

a 



M 
2rc 



(13) 



(14) 



absorbs the singular term analytically, and up is a regular 
correction term that accounts for the effects of the boost. 
In terms oi up the Hamiltonian constraint ((Sj becomes 



V2 



Up 



-I3{l + aupy 



(15) 



CP Aij 
ij ^CP • 



(16) 



Since A^p scales with P, the leading order term in H15(l 
scales with ep and all odd-order terms in ep must vanish. 
An analytical solution to order ep is given by 



Up 



Mel 



8(X + 2rc) 
U2 {rc) f2(cos 



{uq (rc) Po(cos6') 



where 



and 



Po(cos^ 



P2(cos( 



3 2 
- cos I 
2 



(17) 



(18) 



(19) 



are Legendre polynomials and where the radial functions 
uq (rc) and U2 (rc) are 

uq (rc) = M'^ + lOA^'Vc + 40A^V^ + SOA^r^ + 80r^, 

(20) 

and 



M 



U2{rc) = ^(42A^Vc +378A^V?, + I3I6AIV?, 
5rQ V 

2156Al2r^ -I- 1536A1r^ -f 240r^ + 



21M {M + 2rc) \n{ 



M 



M + 2rc- 



(21) 



(see Appendix lA II as well as l29'] ) . Expanding the last 
term of it2 (rc) about rc = shows that it2(rc) = 
0{rh). 

^From the conformal factor (|13l) we can now compute 
several quantities of interest. A short calculation in ap- 
pendix lX4l shows that the ADM energy E of the solution 
is 



5 c p2 

E = M + -Mel + 0{e%)^M + — + 0{ej,). (22) 



We would also like to know the irreducible mass 



M 



A 

16^' 



(23) 



where A is the proper area of the black hole's event hori- 
zon [3^. Since we cannot determine the location of the 
event horizon from initial data alone, we approximate 
this area by the area of the apparent horizon. One might 
expect that computing the irreducible mass to second or- 
der in ep requires the position of the horizon to second 
order as well. As it turns out, however, the second order 
term in the horizon's position cancels out (see Appendix 
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FIG. 1: The top panel shows the ADM energy (solid line) 
and irreducible mass (dashed) for a boosted black hole. The 
lines are drawn according to equations 1221 and 12511 : the 
points are numerical results, with error bars smaller than the 
points. The bottom panel shows the difference between our 
analytical and numerical results. The dotted lines have slopes 
of 4, indicating that the errors in our analytical results are of 
order £% . The deviation from this scaling for very small values 
of P/M is caused by the truncation error in the numerical 
data. 



FIG. 2: The top panel shows perturbative (solid line, 
Eq. H27|l ') and numerical (dots) results for the ADM energy 
of a boosted black hole. The dashed line shows the ADM en- 
ergy of a properly boosted Schwarzschild black hole, Eq. 1281 1. 
The bottom panel shows the deviation (iJnum — Esch)/M in- 
troduced by representing the boosted black hole with con- 
formally flat, maximally sliced puncture data (triangles) as 
well as the error |i?ana — S'numl/M introduced by our pertur- 
bative solution (squares). The dotted lines have slopes of 4, 
indicating that all errors scale as e%. 



lA 3|) . so that we only need the first order correction (see 
Appendix I A 2|) 

h = ^-^epcos9 + 0{el). (24) 

As we show in Appendix lA 31 the irreducible mass M is 
then given by 

M^M + ImeI + 0{e%) =M + ^ + 0{e%). (25) 
8 8M 

We have verified the results for the ADM energy E 122|l 
and irreducible mass M 1)25(1 by comparing them with 
numerical results obtained with a pseudo-spectral ellip- 
tic solver jsll . Spectral techniques are well suited for this 
purpose, since they provide sufficient accuracy to resolve 
the higher-order terms neglected in the analytic treat- 
ment. We note that the spectral elliptic solver appears to 
handle puncture data well, despite initial concerns about 
low differentiablility of u at the punctures. This obser- 
vation would benefit from further investigation. 

In FigHwe show E and M as a function of P, including 
both perturbative and numerical values and their differ- 
ence. Since the leading order error of the perturbative 
results scales with ep we find remarkably good agree- 
ment even for moderately large values of the momentum 



P. For conformally flat Bowen-York based initial data, 
P/M w 0.6 at the predicted location of the innermost 
stable circular orbit (ISCO, e.g. 0, Q), at which point 
the relative error of the perturbative values is approxi- 
mately two percent for the ADM energy, and less than 
one percent for the irreducible mass. There is general 
consensus that the ISCO occurs at a somewhat larger bi- 
nary separation (e.g. |lOl|). where both P and our errors 
are even smaller. 

In (|22l) the ADM energy E is given in terms of the bare 
mass A4 , which has a physical meaning only in the limit 
of infinite separation. It is more intuitive to express E in 
terms of the black hole's irreducible mass. Inverting H25|) 
we have 

M = M-^ + Oie%), (26) 
which we can now insert into (|22|) to find 

E = M+^ + Oi4). (27) 

This result is not surprising but reassuring. As expected, 
we can interpret the ADM energy as the sum of a "rest 
mass" - identified with the irreducible mass - and a New- 
tonian kinetic energy term. Given that we only work to 
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order ej, we expect to find these Newtonian expressions 
only. It is worth emphasizing, however, that these New- 
tonian expressions emerge only after carefully taking into 
account the effect of the boost on the black hole's irre- 
ducible mass, which is an intrinsically relativistic object. 

In Fig. 121 we show perturbative and numerical results 
for the ADM energy as function of P/M . We also include 
the special relativistic result for the ADM energy of a 
properly boosted Schwarzschild black hole, 



ch 



2M 



8M3 



(28) 



Comparing the values for Sana, -Enum, and -Esch we see 
that the error in the ADM energy introduced by con- 
structing transverse-traceless initial data perturbatively 
instead of exactly is only somewhat larger than the devi- 
ation of the ADM energy of conformally flat transverse- 
traceless puncture data from i?sch- Both deviations 
scale as and represent an excess energy over the 
Schwarzschild value that can be interpreted as being as- 
sociated with gravitational radiation. 



B. A Static Black Hole with a Companion 

We would like to derive a result as familiar as (|27|l for a 
static black hole with bare mass i in the presence of a 
second static black hole with bare mass A^2 a coordinate 
distance s away. For this system, the exact solution to 
the Hamiltonian constraint |(SJ) is 



Ml M2 



2rci 2rc2 



(29) 



In equation (|29|) is the coordinate distance from the 
center of hole Ali, and is the coordinate distance 
from the center of hole M.2- Evaluating the ADM energy 
(|A33|) for the conformal factor H29() yields 



perturbation of the horizon and may take into account 
only the zeroth order term 



We then find the irreducible mass 



Ml 



Ml -.^^ del) 



(32) 



(33) 



E = Ml +M2- 



(30) 



(See Appendix IB 2|l . Interchanging indices yields the ir- 
reducible mass of black hole M2- In Figure Owe show 
the irreducible masses of static black holes with compan- 
ions as a function of binary separation s, including both 
perturbative and numerical values and their difference. 
We again find remarkably good agreement to very small 
binary separations. As shown in [20|) the expression H33(l 
turns out to be accurate to at least third order in e^, and 
our comparison in Fig. O demonstrates that the leading 
order error term scales with e^. As computed from the 
transverse-traceless decomposition, for equal masses the 
ISCO occurs at s/Mtot ~ 2.25 (e.g. 4']), where the rel- 
ative error is approximately 10"^. This again sets an 
upper limit, since more realistically the ISCO is believed 
to occur at a somewhat larger binary separation. 

As for the single boosted black holes we again invert 
(123 

M,=M,-^ + 0{e% (34) 

and insert this expression into (|3()|l to find the ADM en- 
ergy E in terms of the irreducible mass 

E^Mi+M2-^^^^ + 0{el). (35) 
s 

This is again the expected Newtonian result which al- 
lows us to interpret the ADM energy as the sum of the 
"rest mass" of the two black holes, identified with the 
irreducible masses, plus the potential energy. 



To calculate the irreducible masses Mi and M2 to the 
order necessary to reproduce the expected classical 
result for E, we approximate i/i in a neighborhood of A^i 
as 



= 1 



Ml 
2r, 



Zs 



(31) 



(see Appendix lB ip . In a neighborhood of A^2, ^ is given 
by the same expression with the indices interchanged. We 
would again like to evaluate the irreducible masses Mi 
and M2 of the two black holes. We proceed exactly as 
we did for isolated boosted black holes, except that in 
this case we only need to expand to order e^. Since the 
correction to the location of the horizon h only enters 
squared into the expression for the irreducible mass (see 
^291) in Appendix I A 3\i . we can completely neglect the 



C. A Binary Black Hole System 

We are now in a position to construct perturbative 
binary black hole solutions by combining the results of 
Sections IIII Al and IIII Bl As it turns out, we can simply 
add both results and obtain a solution that is accurate 
to order 

e = es « ep (36) 

(where we assume that eg and ep are similar for both 
black holes.) More systematically, however, we proceed 
as follows. 

First, we would like to solve the Hamiltonian constraint 
(O for two black holes: one at coordinate location Ci 
with bare mass Mi and momentum Pi, and the other 
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FIG. 3: The top panel shows the irreducible mass for a black 
hole with a companion a coordinate distance s away. Two 
cases with the same total bare mass A^tot = M\-\-M.2 but dif- 
ferent reduced bare masses /ibare = MiMij Mtot are shown; 
M\ = M'2, — 1 and M\ =_§, M-i — \. The line drawn is 
the prediction of equation Il33i , and is the same for each case. 
The points are numerical results: squares for the equal mass 
case, and circles and triangles for the larger and small masses, 
respectively. The bottom panel shows the difference between 
our analytical and numerical results. The dotted lines have 
slopes of —5, indicating that the errors in our analytical re- 
sults are of order t'i. 



at coordinate location C2 with bare mass M.2 and mo- 
mentum P2. As in Sections 1111 Al and 1111 BI we start with 
the corresponding static solution 



Ml M2 



2rci 2rc, 



(37) 



and find a correction u using the puncture method of |22 
for non- vanishing A*-' . Defining 



^ — u + 



1 



with 



Ml M2 



1 

- = 1 , 
a 2rci 2rc2 



(38) 



(39) 



the Hamiltonian constraint Q becomes a regular equa- 
tion for the correction u 



V^u = -/3(1 + au)" 



Here we have again used 



1 



(40) 



(41) 



and Aij is now given by equation (|SJ. Equation H40() 
is nonlinear, and is usually solved numerically. We will 
solve it perturbatively, using our solutions from sections 
IIII Al and 1111 Bl Intuition suggests, and Appendix IC II 
proves that to fourth order in the momenta, u is sim- 
ply the sum of the boost corrections for each hole taken 
separately, 



Up, + Up, +0(6"^), 



(42) 



where up-, and up, are the isolated black hole perturba- 
tions given by equation (flTI . The first neglected terms 
are of order w Cp because, as in Section IIII Al A^^ 
scales with P, so that the leading order term in H40() 
scales with ep and all odd-order terms in ep must van- 
ish. The conformal factor is then 



Ml 



Mo 



2rci 2rc, 



+ up,+up,+0{e'^), (43) 



and it follows immediately that the ADM energy of the 
system to this order is 



cp2 rp2 

E = Mi + -^+M2 + §Th + 0{e^). 



8M1 



8M2 



(44) 



Finding the irreducible mass of each hole again re- 
quires finding the apparent horizons surrounding each 
hole. Since we only consider the leading order perturba- 
tions of a Schwarzschild black hole, we can simply add 
the contributions from the boost and the companion and 
find 



hi 



Ml 
2 



Pi 



16M_ 



■cos 6*1 + 0(5) 



(45) 



where 9i is measured from Pi. We note in particular 
that M2's boost affects Mi only at higher order. The 
irreducible mass of the black hole A^i is then 



Ml = Ml 



M1M2 



8M1 



2s 



+ 0(e') (46) 



(and similar for M2', see Appendix IC 21 1 The correction 
to the irreducible mass of each black hole is the sum of 
the separate corrections for its boost and the presence 
(but not the momentum) of its companion. As before we 
can solve for the bare masses. 



P? M1M2 , o. 
Ml = Ml - ^ + 0{(^) 



8M1 



2s 



(47) 



and similar for M2, and express the ADM energy E in 
terms of the irreducible masses 



Ml + M2 



p2 



M1M2 



+ 0(e'). (48) 



2Mi 2M2 ! 
Again, this result is not surprising but reassuring 
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D. Quasicircular Orbits 

Because of the circularizing effects of gravitational ra- 
diation, binary black holes at reasonably small binary 
separation are expected to be in approximately circu- 
lar orbit. We can construct such systems by minimizing 
the binding energy while keeping the orbital angular mo- 
mentum and irreducible masses fixed (see, e.g., |4fl| | for 
an illustration.) Since Eq. represents the Newtonian 
energy, we will recover expressions of Newtonian orbits in 
what follows, but we also keep track of the error term. We 
define the binding energy as the difference between the 
ADM energy E and the black holes' irreducible masses. 



Eb = E - Ml - M2. 
We also define the total mass 

Mtot = Ml + M2 
and the reduced mass 



M1M2 
Ml + M2 ■ 



(49) 



(50) 



(51) 



In a reference frame where the total momentum is zero 
both individual momenta Pi and P2 must have the same 
magnitude 



Pl=P2 



P. 



(52) 



In such a frame we can then write the binding energy 
(P|l as 



E, = ^-^ + 0{e'). 
2^ s 



(53) 



Not surprisingly, this is again the Newtonian expression 
for the binding energy, and as a consequence the entire 
following discussion is essentially Newtonian. 
Quasicircular orbits are those satisfying 



dEb 
ds 



0, 



(54) 



Mi.M-i,.} 



where J = Ps is the orbital angular momentum (note 
that J = 0(6^1/2) y terms of J, 



Eh 



(55) 



Taking the derivative, we find that quasicircular orbits 
are those for which 



J 



-1/2 



or 



P 



1/2 



+ 0(ei/2) 



+ 0(e3/2) 



(56) 



(57) 
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FIG. 4: The deviation in binding energy iJ(,//x of an equal 
mass binary between our calculation (which yields the Newto- 
nian result) and more sophisticated approaches as a function 
of orbital frequency. We compare with the post-Newtonian 
expansion of |41J as well as numerical relativity results 6, lOj. 
The line for "IPN" has precisely the slope of our error term, 
(Aftotf2)*^^, confirming that our analysis finds the correct 
scaling. 



This condition is equivalent to a Virial relation and justi- 
fies our relation Consequently, quasicircular orbits 
have a binding energy 



E, 



bQC 



IMt, 



2 s 



(58) 



as one might have expected. The orbital angular fre- 
quency, Q, measured at infinity, is 



dEb 
dJ 



Mi,M2,s 

Inserting the binding energy 1)551) we find 



Mtatn = 



(59) 



(60) 



which we recognize as Kepler's third law. Substituting 
into Eq. 158(1 . we finally find 



E, 



bQC 



'2 (Mtot^) 



2/3 



O 



(61) 



Evidently, our expressions for the angular momentum, 
the binding energy and the orbital angular frequency are, 
as expected, simply the Newtonian point-mass expres- 
sions. Consequently, our perturbative approach gives 
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asymptotically correct results at sufficiently large sep- 
arations. The errors in our approach are dominated by 
the first post-Newtonian contributions. As an illustra- 
tion, Fig. ^ shows the difference between our result for 
the binding energy of the binary, Eq. H61|l . and the re- 
sults of more elaborate schemes, namely post-Newtonian 
expansions 1411 and numerically constructed BBH se- 
quences 0, llCl| |. The errors in our expressions scale pre- 
cisely as expected as (Mtot^^)^^^ ~ (Aftot/s)^. 

The graph demonstrates that for a given binary sepa- 
ration the deviations in Ei,/ fi are larger than the devia- 
tions for the individual contributions of the boost and the 
companion, cf. Figs. |2 and |2| This indicates that those 
0(6^) terms which we neglected when constructing the 
binary (cf. Sec. IIII 0|) have larger coefficients than those 
for single boosted or tidally deformed black holes. 



IV. CONSTRUCTING INITIAL DATA 

In this Section we briefly describe how, for arbitrary 
irreducible masses Mi and M2 and binary separation s, 
we can use the results of the previous Sections to con- 
struct approximate binary black hole initial data that 
are accurate to order e. 

Without loss of generality we assume the orbital plane 
to be the 2 = plane of a Cartesian coordinate system, 
and assume the a;-axis to connect the two holes. For 
circular orbits the momenta Pi, 2 must then be aligned 
with the y-axis. In a reference frame where the total 
momentum is zero we must also have 

Pi = -P2. (62) 
^From (|57|l the magnitude of the momenta is given by 



x/s/Mtot 
Using H47|l we find the bare mass 



(63) 



Ml ^Mi{l 



8s (Ml + M2) 2s 



+ 0(e2) (64) 



and is given by the same expression with the indices 
interchanged. To place the center of (irreducible) mass 
at the origin of a Cartesian coordinate system, we choose 



Ci = 



Ml + M2 



for Ml and 



Mis 



Ml + M: 



,0,0 



■,0,0 



(65) 



(66) 



for M2- With the momenta and positions of the two 
black holes given, the extrinsic curvature Aij can now be 
computed from ||SJ) and (T)) . Finally, the conformal factor 
-0 is given H43|) and H17|l . This completes the construction 
of our approximate binary black hole initial data. 



V. SUMMARY 

We construct approximate but analytical binary black 
hole initial data. Adopting the puncture method (see 
[23. I23L l2^ ) we solve the constraint equations of gen- 
eral relativity in the transverse-traceless decomposition 
for perturbations of isolated, static Schwarzschild black 
holes, caused by a boost and the presence of a binary 
companion. We then use a superposition of the indi- 
vidual perturbations to construct analytical black hole 
binary solutions that are accurate to first order in the 
square of the momenta and the inverse of the binary sep- 
aration. In particular, our initial data become exact in 
the limit of infinite separation. 

For binary black hole solutions constructed from 
the transverse-traceless decomposition of the constraint 
equations, quasi-circular orbits are usually identified by 
locating a minimum of the binding energy along se- 
quences of constant mass and angular momentum (see 
Section IIII D|l . For numerical solutions this requires a 
large number of iterations, making the problem quite 
involved. In contrast, for our analytical solutions this 
becomes very simple. In Section IIVI we describe how 
approximate binary black hole solutions in quasicircular 
orbit, for arbitrary mass ratio, can be set up analyti- 
cally. Given the amount of work required to construct 
numerical solutions, our approximate solutions might be 
an attractive alternative for some applications, especially 
for large binary separations where our approximation be- 
comes increasingly accurate. 

Given that we expand only to first order in the square 
of the momentum and the inverse of the binary separa- 
tion, all expressions for global quantities like the energy, 
angular momentum and orbital frequency take their New- 
tonian values. This does not mean, however, that our cal- 
culation is Newtonian. In fact, the above quantities only 
take the expected Newtonian values after carefully tak- 
ing into account corrections to the black holes' irreducible 
mass, which is an intrinsically relativistic quantity. This 
study was partly motivated by our curiousity in the be- 
havior of the bare mass along sequences of constant irre- 
ducible mass (compare footnote [35] and the surrounding 
discussion in M5|.) Moreover, our solutions provide ap- 
proximate values for the conformal factor and the extrin- 
sic curvature together with the above global quantities, 
thereby specifying complete sets of initial data describing 
binary black holes in approximate quasicircular orbit. 
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APPENDIX A: A BOOSTED BLACK HOLE 

1. The Conformal Factor 

We would like to solve the Hamiltonian constraint Q 
for a boosted black hole with bare mass A4 and momen- 
tum P located at C. We follow the puncture approach 
of and write the conformal factor ip as 



1 



Up 



(Al) 



(see Section lill A|) . in which case the problem reduces to 
solving 



(A2) 



for the correction term up. Here the coefficients a and 
(3 are 



a = 1 



M 

2rc 



(A3) 



and 



/3=-aMJ^-A^p. (A4) 
The extrinsic curvature Af^^ is given by Q and yields 

(A5) 



aCP 

^ij ^CP 



9P2 fi 



' c 



cos^ 



where 9 is measured from P. Since P enters into the 
equations only squared, all odd order terms in P must 
vanish. We construct a perturbative solution that is ac- 
curate to order P^, so that the leading order error term 
scales with P'^. Up to this order, (|A2|) reduces to 



V'^up+(3 + 0{e%) = 



(A6) 



where €p — Pj M.. Because of the axisymmetry of this 
problem, it is reasonable to guess a solution of the form 

up = e^p ({2o(rc)Po(cos^) + U2(rc)P2(cos ^)) + 0(e|,), 

(A7) 

where Po(cos6') = 1 and P2(cos6') = (3cos2 6' - l)/2 are 
Legendre polynomials. Solutions satisfying the boundary 
conditions up ^ as re ^ oo and (9wp/9rc)|rc=o ~ 
are 



Wo (»'c) 



M 



8(yW + 2rc) 



-A 



(A8) 



and 



1316A^V^ + 2156M^r^ + 1536A^r^ 



240rp 



21M {M + 2rcf\n{ 



M 



M + 2rc- 



(A9) 



(see also [23). In Section [ill Al we have factored out 
terms that are common to both uo(^c) and ^2(^0) (see 



2. The Apparent Horizon 

In axisymmetry, the apparent horizon is located at a 
distance rc{0,(j)) = h{9) from a point C located on the 
axis of symmetry, where hiff) satisfies the following ordi- 
nary differential equation |43l | 



dedgh = -ricmAU^'u^ -i^y 

an 

, ds , 



subject to the boundary condition 
dh. 



de 



0=0, TT 



0. 



In (|IlO)l we use the vectors 

rrii = (1, -deh, 0), 

and 

u' = {deh, 1,0). 

We also abbreviate 

,ds 



2 



(AlO) 

(All) 

(A12) 
(A13) 

(A14) 



where capital letters A, B, etc. run over the coordinates 
re and 0, but not (j). All coefficients are evaluated at 
the horizon location, rc = h. Given our assumption of 
conformal flatness we also have 



1 



r^, sin^ 



(A15) 



The black hole's momentum enters (jA10|l through the 
extrinsic curvature Kij. 

It is easy to verify that for vanishing momentum P — Q 
the well-known isotropic horizon radius 



M 
2 



(A16) 
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satisfies (|A10|) and the boundary conditions. To find a 
perturbative solution for non-zero P we expand 



h ^ ho + hp + hp2 + 0(ep), 



(A17) 



where hp and hp2 are corrections of first and second 
order in ep, respectively. 

To find the first order correction h p we need to express 
all terms in IjAlOfl up to order cp, which yields 

dedeh dedehp + 0{el), 

(^)V^r^,m^ = icot0)dehp-^ + O{el), 

il^'^r'^'fyu^u^KAB = ^-'Kgg + 0{el) (Al8) 
cLu 



3P 
'32 



cos 61 + 0(61,) 



, ds . 



= -^cos0 + O(e|,). 

The resulting equation for hp is 

3 p 

d9dehp~hp+(cotd)dehp~—cos0+O(el) =0. (A19) 

16 

The solution satisfying the boundary conditions is 



hp = cos6' + 0(e?,) 
16 



(A20) 



(compare |23|). One might expect that in order to find 
the irreducible mass to order ej, we also need the next 
order term hp2. However, as we will see in Section lA 31 
this second order terms cancels out, so that in fact we 
only need the first order correction hp. 



3. The Irreducible Mass 

We would now like to find the irreducible mass of a 
boosted black hole, to second order in ep. We approxi- 
mate the irreducible mass as 



M = 



A 

16^' 



(A21) 



where A is the proper area of the apparent horizon of the 
black hole js^. It can be shown (see, e.g.. Appendix D 
in lill) that 



A: 



(A22) 

where -0 and rc are evaluated at the horizon location, 
rc — h. From the expansion IjAlTp of h we have 




dh\^ / dhp^ ^ 



de \ 89 



(A23) 



since ho is spherically symmetric. The horizon area HA22|) 
therefore splits into the two terms 



A = I I sin 9d9d(t> + 



(A24) 



^0 

2n pTT 



,4 f 9hp 



/o JO V 99 
We now insert the expansions 
M 



sin 9d9d4) + 0{el). 



and 



ij = l + — +up + Oie%) 
Ire 



M 



(A25) 



rc = h = '-^ + hp + hp2 + 0{e%). 



(A26) 

Taking these terms to their respective powers we find 



V' = 16 1 - 



4hp 4hp2 14h^ 



^ + 2up ) + 0(e|,), 
(A27) 



and 



rl ^ — +Mhp + Mhp2 +h% + 0{el), (A28) 



which we can now insert into (|A24|) to find 

A = WttM'^ + WttM'^ I upsm9d9 
Jo 

-flGTT^ hlsm9d9 + Tr ^4 | sin9de 



(A29) 



Quite remarkably, the second order term hp2 has can- 
celed out of this expression. 

To evaluate the horizon area we insert U p (eq uation 
(jA?!) ) fro m Ap pendix I A H and hp (equation (A20P '1 from 
Appendix IA 21 Since we are only working to second order 
in ep it is sufficient to evaluate up at the unperturbed 
horizon location /iq: 



p2 



Up 



:(- 3422-|- 6721n2- 



25607W2 

(11196 - 2016 In 2) cos^ 61) + 0(e|,). (A30) 



In yet another remarkable happenstance the unattractive 
log terms disappear when the integration is carried out, 
and we are left with 

A=167:M\l + ^) + 0{e%), (A31) 



and therefore 



M = M{l + ^) + 0{e%). (A32) 

Even though we have carried out this calculation only to 
order ep, it is clear that neither the area nor the mass 
can depend on the sign of P, so that all odd order terms 
must disappear. We have also verified this by comparing 
with numerical data in Fig. ^ 
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4. The ADM Energy 



The Irreducible Mass 



The ADM energy E is defined as 



E 



1 

'2^ 



(A33) 



According to (|13() we write the conformal factor as '0 = 
1 + M/2rc + Up. The correction term up is given by 
(IA7(I . and evaluating its leading order term as re — > c)o 
we find 



5P^ 



Up 



IGMrc 

The ADM energy therefore becomes 

c p2 



(A34) 



(A35) 



APPENDIX B: A STATIC BLACK HOLE WITH A 
COMPANION 

1. The Conformal Factor 

The exact solution for the conformal factor describ- 
ing two static black holes is given by 



1 



Ml M2 



(Bl) 



Focusing on Ml we would like to eliminate the depen- 
dence on , which we can do in terms of the expansion 



1 

rc2 



^ 00 



s — ' \ s 

n=0 



P„ f cos ( 



(B2) 



Here 9 is the angle between rci and T12 = C2 — Ci. Since 
we are only interested in terms up to order Cs = M / s it 
is sufficient to keep the n = term as long as we restrict 
analysis to a neighborhood oi Mi, where rci ~ Mi- 
Near A^i we then have 



1 1 



rc2 



and hence 



2rci 2s 



(B3) 



(B4) 



and similar for M2- Borrowing the notation for boosted 
black holes we denote 



Zs 



(B5) 



Evidently the effect of a companion on an otherwise 
spherical black hole must be axisymmetric in nature. To 
find the irreducible mass of the perturbed black hole M 1 
we can therefore again evaluate the integral ljA29p . with 
Up replaced by uci (equation (|B5|I .) As before we can 
expand the horizon location as 



h ^ ho + hs + hs2 + O(e^). 



(B6) 



Unlike in the case of a boosted black hole, however, we 
are only interested in terms up to order e^. Since the 
leading order term scales with and enters into the 
irreducible mass only squared, we may neglect its con- 
tributions. The only terms that remain in (|A29|I are 
therefore the first two. Inserting (|B5|I we immediately 
find 



and hence 



Ai = l6TTMlil + —) + 0{el), 



Mi^Mi + ^ + Oi.l). 
Zs 



(B7) 



(B8) 



Interchanging indices yields M2- 

APPENDIX C: A BINARY BLACK HOLE 
SYSTEM 

1. Conformal factor 

We would now like to solve the equation 

V^u = -/3 + ©(e^) (CI) 

for two boosted black holes (here e is defined in H36|) ). 
The source term P is 



(C2) 



where a is defined in H39|) and the extrinsic curvature 
Aij is given by © as the sum of the individual terms for 
the two black holes. Given that this source term is non- 
linear, it is not immediately evident that we can simply 
superpose two single boosted black hole solutions 



u = + 7iP2 + 0{e^), 



(C3) 



each of which satisfies (IC1|) with their individual source 
term, e.g. 



(C4) 



The superposition IjCSp only satisfies (jCip to the required 
order if 



and similar for M2- 



p-Pi-f32^0{e'). (C5) 
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This, however, is indeed the case. Any point in space is 
at least a coordinate distance s/2 separated from one of 
the black holes. For concreteness, consider a point that 
is at least a distance s/2 separated from A42- Expanding 
/? at this point shows 

P = (3i+ 0{e^), (C6) 

and we also have f32 = 0{e^). Therefore, the superposi- 
tion l|C3l) is a solution to ljCl|) to the required order. 

We point out that the non-linearity of /3 prevents us 
from also including spin in our analysis. Allowing for 
non-zero spin amounts to including additional terms in 
the extrinsic curvature, which then lead to cross-terms 
in the source term (|C2|I . Perturbative solutions describ- 
ing stationary, spinning black holes exist, but because 
of these cross-terms they cannot simply be added to the 
perturbative solutions describing non-spinning, boosted 
black holes to construct spinning, boosted black holes. 

2. Irreducible mass 

We finally argue why the correction to the irreducible 
mass of a boosted black hole with a companion is the 
sum of the corrections for a boost and the presence of 



a companion. In the neighborhood oi Mi, for example, 
the conformal factor is given by H43|l . 

ij = l + ^^+uc,+up„ (C7) 

where = A^2/(2s) + C(e^), and MiS horizon is lo- 
cated at 

rc,^hi = ^ + hp + Oie). (C8) 

As before, higher order corrections to the location of the 
horizon drop out. We can now proceed exactly as in 
Appendix I A 31 using the combined correction v = uc\ + 
up-^ to the conformal factor, and find 

Ai = 167tMI{1 + ^ + + 0{e'). (C9) 

The irreducible mass is therefore 

^ -^^(^ + 8^4? + + ^^''^ ^^^^^ 
as expected. 
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